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Abstract 

We investigate experimentally wind tunnel turbulence generated by multiscale/fractal grids pertaining 
fN| | to the same class of low-blockage space-filling fractal square grids. These grids are not active and 

$_i ' nevertheless produce very much higher turbulence intensities u'/U and Reynolds numbers Re\ than 

higher blockage regular grids. Our hot wire anemometry confirms the existence of a protracted pro- 
duction region where turbulence intensity grows followed by a decay region where it decreases, as first 
reported by Hurst & Vassilicos (l5| . We introduce the wake-interaction length-scale x* and show that 
the peak of turbulence intensity demarcating these two regions along the centreline is positioned at 
about 0.5x*. The streamwise evolutions on the centreline of the streamwise mean flow and of various 
statistics of the streamwise fluctuating velocity all scale with x*. Mean flow and turbulence intensity 
profiles are inhomogeneous at streamwise distances from the fractal grid smaller than 0.5x*, but appear 
quite homogeneous beyond 0.5:e*. The velocity fluctuations are highly non-gaussian in the production 
region but approximately gaussian in the decay region. Our results confirm the finding of Seoud & 
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Vassilicos [31| that the ratio of the integral length-scale L u to the Taylor microscale A remains con- 



stant even though the Reynolds number Re\ decreases during turbulence decay in the region beyond 
0.5x*. As a result the scaling L u /X ~ Re\ which follows from the u' 3 /L u scaling of the dissipation 
rate in boundary-free shear flows and in usual grid-generated turbulence does not hold here. This 
extraordinary decoupling is consistent with a non-cascading and instead self-preserving single-length 



, scale type of decaying homogeneous turbulence proposed by George & Wang [121], but we also show 

that L u /X is nevertheless an increasing function of the inlet Reynolds number Re®. Finally, we offer a 
detailed comparison of the main assumption and consequences of the George & Wang theory against 
QO . our fractal-generated turbulence data. 
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G\ ; 1 Introduction 
O 

Which turbulence properties are our current best candidates for universality or, at least, for the defi- 
nition of universality classes? The assumed independence of the turbulence kinetic energy dissipation 
rate on Reynolds number Re in the high Re limit is a cornerstone assumption on which Kolmogorov's 
phenomenolo gy i s built and on which one-point and two-point closures and LES rely, whether directly 
or indirectly [3], jjj, (29| . |3oT ] . 21 1. This cornerstone assumption is believed to hold universally (at 



least for relatively weakly strained/sheared turbulent flows). It is also related to the universal tendency 
of turbulent flows to develop sharp velocity gradients within them and to the apparently universal 



geometrical statistics of these gradients [32j, as to the apparently universal mix of vortex stretching 
and compression (described in some detail by Tsinober [35] who introduced the expression "qualitative 
universality" to describe such ubiquitous qualitative properties). 

Evidence against universality has been reported since the 1970s, if not earlier, in works led by 



Roshko, Lykoudis, Wygnanski, Champagne and George (see for example [111 ] and references therein as 



well as the landmark work of Bevilaqua & Lykoudis 0] and more recent works such as (3r3 | and |2 
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to cite but a few) and has often been accounted for by the presence or absence of long-lived coherent 
structures. Coherent/persistent flow structures can actually appear at all scales and can be the carrier 
of long-range memory, thus implying long-range effects of boundary/inlet conditions. 

In summary, kinetic energy dissipation, vortex stretching and compression, geometrical alignments, 
coherent structures and the universality or non- universality of each one of these properties are central to 
turbulent flows with an impact which includes engineering turbulence modelling and basic Kolmogorov 
phenomenology and scalings. Is it possible to tamper with these properties by systematic modifications 
of a flow's boundary and/or inlet/upstream conditions? 

To investigate such questions, new classes of turbulent flows have recently been proposed which 
allow for systematic and well-controlled changes in multiscale boundary and/or upstream conditions. 
These new classes of flows fall under the general banner of "fractal-generated turbulence" or "multiscale- 
generated turbulence" (the term "fractal" is to be understood here in the broadest sense of a geometrical 
structure which cannot be described by any non-multiscale way, which is why we refer to fractal and 
multiscale grids interchangeably). These flows have such unusual turbulence properties (l5| . (3l| that 
they may directly serve as new flow concepts for new industrial flow solutions, for example conceptually 
new energy-efficient industrial mixers [6(. These same turbulent flow concepts in conjunction with 
conventional flows such as turbulent jets and regular grid turbulence have also been used recently 
for fundamental research into what determines the dissipation rate of turbulent flows and even to 
demonstrate the possibility of renormalising the dissipation constant so as to make it universal at 
finite, not only asymptotically infinite, Reynolds numbers (see 25 1, [l3|). These works have shown that 
the dissipation rate constant depends on small-scale intermittency, on dissipation range broadening 
and on the large-scale internal stagnation point structure which itself depends on boundary and/or 
upstream conditions. In the case of at least one class of multiscale-generated homogeneous turbulence, 
small-scale intermittency does not increase with Reynolds number [33j and the dissipation constant is 
inversely proportional to turbulence intensity even though the energy spectrum is broad with a clear 
power-law shaped intermediate range [31| . [15l |. In this paper we investigate this particular class of 
multiscale-generated turbulent flows: turbulent flows generated by low-blockage space-filling fractal 
square grids. 

Grid-generated wind tunnel turbulence has been extensively investigated over more than seventy 
years [3| and is widely used to create turbulence under well controlled conditions. This flow has the great 
advantage of being nearly homogeneous and isotropic downstream However, its Reynolds number 
is not large enough for conclusive fundamental studies and industrial mixing applications. Several 
attempts have been made to modify the grid so as to increase the Reynolds number whilst keeping 
as g ood homogeneity and isotropy as possible: for example jet-grids by Mathieu's (2^ and Corrsin's 
10 1 groups (who may have been inspired by Betchov's porcupine non-stationary, so-called active, 



grids by Makita [22| followed by Warhaft's group [27| and others, and most recently passive grids with 
tethered spheres attached at each mesh corner [38|. Jet-grids and active grids have been very successful 
in increasing both the integral length-scale and the turbulence intensity whilst keeping a good level 
of homogeneity and isotropy. The three different families of fractal/multiscale grids introduced by 
Hurst & Vassilicos [15l | generate turbulence which becomes approximately homogeneous and isotropic 
considerably further downstream than jet-grids and active grids, but achieve comparably high Reynolds 
numbers even though, unlike jet-grids and active grids, they are passive. However, the most important 
reason for studying fractal/multiscale-generated turbulence is that it can have properties which are 
clearly qualitatively different from properties which are believed to be universal to all other grid- 
generated turbulent flows and even boundary- free shear flows for that matter. 

In this paper we report the results of an experimental investigation of turbulent flows generated by 
four low-blockage space-filling fractal square grids. The grids used in our study are described in the 
next section and the experimental set up (wind tunnels and anemometry) is presented in section 3. Our 
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results are reported in section 4. Specifically, in subsection 4.1 we introduce the wake-interaction length- 
scale and use it to derive and explain the scaling of the downstream peak in turbulence intensity 
which was first reported in (l5| . We also show in this subsection that the streamwise dependence of 
the streamwise turbulence intensity is independent of inlet velocity and fractal grid parameters if x* 
is used to scale streamwise distance. In subsection 4.2 we confirm the far field statistical homogeneity 
first reported in 31| and, for the first time, present near-field profiles illustrating the evolution from 
near-field inhomogeneity to far-field homogeneity. Subsection 4.3 contains a detailed report on the 
skewness and flatness of the fluctuating velocities illustrating how they become gaussian in the far field 
following a clearly non-gaussian near-field behaviour which peaks at 0.2x*. Finally, in subsection 4.4 
we report a significant improvement and generalisation of the single-scale self-preservation theory of 
George & Wang fl~2| which shows that there are many more single-scale solutions to the spectral energy 
equation than originally thought. Subsections 4.5 to 4.8 make use of this multiplicity of solutions for an 
analysis of our data that is significantly finer than in previous studies of fractal-generated turbulence 
and which confirms the self-preserving single-scale nature of the far-field decaying fractal-generated 
turbulence in terms of the behaviours of the integral scale, the Taylor microscale, the energy spectrum 
and the turbulence intensity. 

Finally, in section 5 we conclude and discuss some of the issues raised by our investigation. 



2 The space-filling fractal square grids 



Turbulent flows are generated in this study by the planar and space-filling multiscale/fractal square 
grids first introduced and described in [15| . The main characteristics of those grids are summarised 
as follows. In general, multiscale/fractal grids consist of a multiscale collection of obstacles/openings 
which may be all based on a single specific pattern that is repeated in increasingly numerous copies at 
smaller scales. For the present work, the pattern used is a empty square framed by four rectangular bars 
as shown in figure 1(a) Each scale- iteration j is characterized by a length Lj and a thickness tj of these 
bars. At iteration j there are four times more square patterns that at iteration j — 1 (1 < j < N where 
N is the total number of scales) and their dimensions are related by Lj = R^Lj-i and tj = R^tj-i. The 
scaling factors Rl and Rt are independent of j and are smaller or equal than 1/2 and 1 respectively. As 
explained in (l5| . the fractal square grid is space filling when its fractal dimension takes the maximum 
value 2, which is the case when Rl = 1/2. 
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Figure 1: (a) Space- filling multiscale/fractal square grid generating pattern. Examples of planar mul- 
tiscale/fractal square grids used in the present work with N=4 scales (b) and N=5 scales (c). 



A total of four different planar space-filling fractal square grids have been used in the wind tunnel 
experiments reported here. The complete planar geometry of these grids is detailed in table [TJ Scaled- 
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down diagrams of two of these grids are displayed in figures 1 (b) and 1 (c) Multiscale /fractal grids are 
clearly designed to generate turbulence by directly exciting a wide range of fluctuation length-scales 
in the flow rather than by relying on the non-linear cascade mechanism for multiscale excitation. The 
latter approach is the classical one and is exemplified by the use of regular grids as homogeneous 
turbulence generators. 



Grid 


SFG8 


SFG13 


SFG17 


BFG17 


Lo(mm) 


237.5 


237.7 


237.8 


471.2 


Li(mm) 


118.8 


118.9 


118.9 


235.6 


Li2{mm) 


59.4 


59.4 


59.5 


117.8 


L^(mm) 


29.7 


29.7 


29.7 


58.9 


L^(mm) 








29.5 


to(mm) 


14.2 


17.2 


19.2 


23.8 


ti(mm) 


6.9 


7.3 


7.5 


11.7 


t2(mm) 


3.4 


3.1 


2.9 


5.8 


£3 (mm) 


1.7 


1.3 


1.1 


2.8 


£4 (mm) 








1.4 



Table 1: Geometry of the space- filling fractal square grids. 



As explained in [151 ] . the complete design of space-filling grids requires a total of four independent 
parameters such as : 

• N : the number of scales (N — 1 being the number of scale- iterations) , 

• Lq : the biggest bar length of the grid 

• to : the biggest bar thickness of the grid 

• ijv_i : the smallest bar thickness of the grid 

The smallest bar length L^^\ on the grid is determined by Rl = 1/2 and N. Note, also, that the 
fractal grids are manufactured from an acrylic plate with a constant thickness (5mm) in the direction 
of the mean flow. 

Hurst & Vassilicos (l5| introduced the thickness ratio t r = - and the effective mesh size 

AT 2 , 

M eff = -jr^^ (1) 

where P is the fractal perimeter length of the grid, T 2 is the tunnel's square cross section and a is the 
blockage ratio of the grid defined as the ratio of the area A covered by the grid to T 2 : 

a =T 2 = 1 Y 2 ' () 

These quantities are derivable from the few independent geometrical parameters chosen to uniquely 
define the grids. When applied to a regular grid, this definition of M e ff returns its mesh size. When 
applied to a multiscale grid where bar sizes and local blockage are inhomogeneously distributed across 
the grids, it returns an average mesh size which was shown in fl~5| to be fluid mechanically relevant. 

A total of four space-filling fractal square grids have been used in the present work. They all have 
the same blockage ratio a = 0.25 (low compared to regular grids where, typically, a is about 0.35 or 0.4 
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or so [8j, [7J]) and turn out to have values of M e ff which are all very close to 26.4mm. Three of these 
grids, refered to as SFG8, SFG13 and SFG17, differ by only one parameter, t r , and as a consequence 
also by the values of tj (0 < j < N — 1) as t r was chosen to be one of the four all-defining parameters 
along with the fixed parameters N = 4, L3 = 29.7mm and a = 0.25. The fourth grid, BFG17, has one 
extra iteration, i.e. N = 5 instead of N = 4 but effectively the same smallest length, i.e. L4 = 29.5mm, 
and a value of t r very close to that of SFG17. It is effectively very similar to SFG17 but with one extra 
fractal iteration. The main characteristics of these grids are summarized in tables 1 and [2] which also 
includes values for L r = Lq/L^-i- 



Grid 


N 


Lo(mm) 


to (mm) 


L f 




a 


M e ff(mm) 


SFG8 


4 


237.5 


14.2 


8 


8.5 


0.25 


26.4 


SFG13 


4 


237.7 


17.2 


8 


13.0 


0.25 


26.3 


SFG17 


4 


237.8 


19.2 


8 


17.0 


0.25 


26.2 


BFG17 


5 


471.2 


23.8 


16 


17.0 


0.25 


26.6 



Table 2: Main characteristics of the fractal square grids. 



In addition to the fractal grids, we have also performed a comparative study of turbulence generated 
by a regular grid, refered as SRG hereafter, made of a bi-plane square rod array. Table [3] presents the 
main properties of this grid. Its blockage ratio is higher than that of our fractal grids and closer to 
the usual values found in literature for regular grids (see e.g. 0|). The regular grid SRG also has 
slightly higher mesh size. 



Grid 


N 


Lo(mm) 


to (mm) L r t r a M e ff(mm) 


SRG 


4 


460 


6 8 1 0.34 32 



Table 3: Main characteristics of the regular grid. 



3 The experimental set-up 
3.1 The wind tunnels 

Measurements are performed in two air wind tunnels, one which is open-circuit with a 5m long and 
T = 0.46m wide square test section and one which is recirculating with a 5.4m long and T = 0.91m 
wide square test section. A generic sketch of a tunnel's square test section is given in figure [2] for the 
purpose of defining spatial coordinate notation. The arrow in this figure indicates the direction of the 
mean flow and of the inlet velocity [Too- The turbulence-generating grids are placed at the inlet of the 
test section. 

The fractal grids SFG8, SFG13, SFG17 and the regular grid SRG were tested in the open circuit 
tunnel whereas the fractal grid BFG17 was tested in the recirculating tunnel. 

The maximum flow velocity without a grid or any other obstruction is 33m/ s in the T = 0.46m 
open-circuit tunnel. Turbulence-generating grids were tested with three values of the inlet velocity [Too 
in this tunnel: 5.2m/s, lOm/s and 15m/s. The uniformity of the inlet velocity at the convergent's 
outlet, checked with Pitot tube measurements, is better than 5%. The residual turbulence intensity in 
the absence of a turbulent-generating grid is about 0.4% along the axis of the tunnel. 

In the T = 0.91m recirculating tunnel, the maximum flow velocity without a grid or any other 
obstruction is 45m/s. The inlet velocity Uoo was fixed at 5.2m/s in this facility when testing the 
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turbulence generated by the BFG17 grid. The entrance flow uniformity is better than 2% and a very 
low residual turbulence intensity (< 0.05%) remains in the test section in the absence of a turbulence- 
generating grid or obstacle. 

In both tunnels, the temperature is monitored during measurement campaigns thanks to a ther- 
mometer sensor located at the end of the test section. The inlet velocity Uoo is imposed by measuring 
the pressure difference in the tunnel's contractions with a micromanometer Furness Controls MCD1001. 



3.2 Velocity measurements 

A single hot-wire, running in constant-temperature mode, was used to measure the longitudinal velocity 
component. The DANTEC 55P01 single probe was driven by a DISA 55M10 anemometer and the probe 
was mounted on an aluminum frame allowing 3D displacements in space. A systematic calibration of 
the probe was performed at the beginning and at the end of each measurement campaign and the 
temperature was monitored for thermal compensation. The sensing part of the wire (PT-O.IRd) was 
5/j.m in diameter (d w ) and about 1mm in length l w so that the aspect ratio was about 200. Our 
spatial resolution l w /rj ranges between 2 and 9 for all the measurements. The estimated frequency 
response of this anemometry system is 1.5 to 9 times higher than the Kolmogorov frequency f v = 
The spatially-varying longitudinal velocity component u(x) in the direction of the mean flow was 
recovered from the time-varying velocity u(t) measured with the hot-wire probe by means of local 
Taylor's hypothesis as defined in fl6| . 

The signal coming from the anemometer was compensated and amplified with a DISA 55D26 signal 
conditioner to enhance the signal to noise ratio which is typically of the order of 45cLB for all measure- 
ments. The uncertainties on the estimation of the dissipation rate e due to electronic noise occurring 
at high frequencies (wavenumbers) is lower than 4% for all our measurements. The conditioned signal 
was low-pass filtered to avoid aliasing and then sampled by a 16 bits National Instruments NI9215 USB 
card. The sampling frequency was adjusted to be slightly higher than twice the cut-off frequency. The 
sampled signal was then stored on the hard-drive of a PC. The signal acquisition was controlled with 
the commercial software LabVIEW©, while the post-processing was carried out with the commercial 
software MATLAB©. 

The range of Reynolds numbers and lengths-scales of the turbulent flows generated in both tunnels 
by all our grids are summarized in table 0] The longitudinal integral length-scale L u was obtained 




Figure 2: Wind-tunnel sketch and coordinate notation 
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by integrating the autocorrelation function of the fluctuating velocity component u(x) (obtained by 
subtracting the average value of u{x) from u(x)): 



(u(x)u(x + A)) 



dA 



(3) 



i (u{x) 2 } 

where the averages are taken over time, i.e. over x in this equation's notation, where x is obtained 
from time t by means of the local Taylor hypothesis. In this paper we use the notation v! = W (u(x) 2 }. 
The Taylor microscale A was computed via the following expression: 



<u*> 



^ {(du/dx)' 



(4) 



Finally, using the kinematic viscosity v of the fluid (here air at ambient temperature) we also 
calculate the length-scale 



1/4 



IhUdu/dxY 



(5) 



which is often refered to as Kolmogorov microscale. We have estimated from the turbulent kinetic energy 
budget that the uncertainties in the computation of A and r\ are lower than 10% and 5% respectively 
for all our measurements. 



T[m] 


Uoo[m/s] 


L u [mm] 


A [mm] 


rj[mm] 


Rex 


0.46 


5/10/15 


43-52 


4-7.4 


0.11 - 0.32 


140 - 370 


0.91 


5 


50-70 


7-10 


0.3 - 0.45 


60 - 220 



Table 4: Main flow characteristics: T is the wind tunnel width, L u is the longitudinal integral length- 
scale, A the Taylor-microscale, rj the Kolmogorov scale and Re\ the Taylor based Reynolds number. 



4 Results 



Hurst & Vassilicos [15j found that the streamwise and spanwise turbulence velocity fluctuations gener- 
ated by the space-filling fractal square grids used here increase in intensity along x on the centreline till 
they reach a point x = x pea k beyond which they decay. Thus they defined the production region as being 
the region where x < x pea k and the decay region as being the region where x > x pea ^. They also found 
that various turbulence statistics collapse when plotted as functions of x/x pea k and they attempted to 
give an empirical formula for x pea ^ as a function of the geometric parameters of the fractal grid. It 
was also clear in their results that the turbulent intensities depend very sensitively on parameters of 
the fractal grids even at constant blockage ratio, thus generating much higher turbulence intensities 
than regular grids. An understanding and determination of how x pea k and turbulence intensities de- 
pend on fractal grid geometry matters critically both for achieving a fundamental understanding of 
multiscale-generated turbulence and for potential applications such as in mixing and combustion. In 
such applications, it is advantageous to generate desired high levels of turbulence intensities at flexibly 
targeted downstream positions with as low blockage ratios and, consequently, pressure drop and power 



input, as possible. An important question left open, for example, in [15j is whether x pea k does or does 
not depend on U^. 
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This section is subdivided in eight subsections. In the first we study the steamwise profiles of the 
streamwise mean velocity and turbulence intensity and, in particular, determine x pea k- In the second 
we offer data which describe how homogeneity of mean flow and turbulence intensities is achieved when 
passing from the production to the decay region. In the third subsection we present results on the 
turbulent velocity skewness and flatness. The fourth and eightth subsections are a careful application 



of the theory of George & Wang [12l | to our data and the fifth, sixth and seventh subsections are an 
investigation of the single-length scale assumption of this theory and its consequences, in particular the 
extraordinary property first reported in (3l| that the ratio of the integral to the Taylor length-scales is 
independent of Re\ = u'X/u in the fractal-generated homogeneous decaying turbulence beyond x pea k- 

4.1 The wake-interaction length-scale x+ 

The dimensionless centreline mean velocity Uc/Uoo and the centreline turbulence intensity u' c /Uc are 
plotted in figures 3(a) and |3(b)| for all space-filling fractal square grids as well as for the regular grid 



SRG. For the latter, we have fitted the turbulence intensity u' c /Uc with the well-known power-law 
A { ^m*® J where the dimensional parameter A, the exponent n and the virtual origin xo have been 
empirically determined following the procedure introduced by Mohamed Sz LaRue (2rj | . Our results are 
in excellent agreement with similar results reported in the literature for regular grids, e.g. n = 1.41 is 



very close to the usually reported empirical exponent (see [261 ] and references therein). 

Figure |3(b)| confirms that a protracted production region exists in the lee of space-filling fractal 
square grids, that it extends over a distance which depends on the thickness ratio t r and that it is 
followed by a region (the decay region first identified in [15]) where the turbulence decays. The existence 



of a distance x pea k where the turbulence intensity peaks is clear in this figure. Figure 3(a) shows that 
the production region where the turbulence increases is accompanied by significant longitudinal mean 
velocity gradients which progressively decrease in amplitude till about after x = x pea k where they more 
or less vanish and the turbulence intensity decays. 

Our data show that the centreline mean velocity is quite high compared to Uoo on the close down- 
stream side of our fractal square grids and remains so over a distance which depends on fractal grid 
geometry before decreasing towards Uoo further downstream. This centreline jet-like behavior seems to 
result from the relatively high opening at the grid's centre where blockage ratio, which is inhomoge- 
neously distributed on the grid, seems to be locally small compared to the rest of the grid. The initial 
plateau is therefore characterized by a significant velocity excess (Uc/Uoo > 1-35). One can also see 



in figure 3(a) that the mean velocity remains larger than Uoo even far away from the grids. We have 
checked that this effect is consistent with the small downstream growth of the boundary layer on the 
tunnel's walls. 

The space-filling fractal grids SFG8, SFG13 and SFG17 are identical in all but one parameter: 



the thickness ratio t r . It is therefore clear from figures 3(a) and 3(b) that t r plays an important role 
because, even though the streamwise profiles of Uc/Uoo and u' c /Uc are of identical shape for SFG8, 
SFG13 and SFG17, Uc/Uoo decreases, u' c /Uc increases and x pea k decreases when increasing t r whilst 
keeping all other independent parameters of the grid constant. 

However, the parameter t r cannot account alone for the differences between the SFG17 and BFG11 
grids. These two grids have the same blockage ratio a and very close values of t r and effective mesh size 
M e ff. What they do mainly differ by, are their values of Lq (by a factor of 2), the number of fractal 



iterations (N = 4 for SFG17, N = 5 for BFG17) and the largest thickness to. Figures 3(a) and 3(b) 
show clearly that when t r is kept roughly constant and other grid parameters are varied (such as Lq), 
then x pea k an d the overall streamwise profiles of Uc/Uoo and of u' c /Uc change in ways not accounted 
for by the changes between SFG8,SFG13 and SFG17. 

Comparing data obtained downstream from different space-filling fractal square grids, Hurst & 
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Figure 3: Streamwise evolution of (a) the centreline mean velocity Uc/Uoo and of (b) the turbulence 
intensity u' c /Uc for all fractal grids with respect to the streamwise distance x scaled by M e ff. Stream- 
wise evolution of (c) the centreline mean velocity Uc/Uoq and of (d) the turbulence intensity u' c /Uc 
normalised by its value at x = x pea k for all fractal grids with respect to the streamwise distance x scaled 



by Xhv = 75t m i n T / L m i n [15j ({Too = 5.2m/ s). For reference, data obtained with the regular grid SRG 

( _ \ -1-41 

and Uoo = Wm/s, are also reported (* symbols in (b)). The power law (u' c /Uc) = A ( X M x ® 
which is also shown (solid line in (b), A = 0.129, xq = 0) fits this SRG data very well. 
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Vassilicos [15l | suggested that the streamwise evolution of turbulence intensity, i.e. x pea k, can be scaled 
by the length-scale xhv = 75 Their empirical formula might appear to account for the difference 

l-'min 

between the SFG17 and BFG17 grids in figure 3(b) because T is double and t m in is larger by a factor 



1.3 for BFG17 compared to SFG17. However, Hurst &; Vassilicos [15j did not attempt to collapse data 
from different wind tunnels, and we now show how such a careful collapse exercise involving both the 
mean flow and the turbulence intensity leads to a different formula for x pea k- 

In figures 3(c) and 3(d) we plot the streamwise evolution of Uc/Uoo and of u'JUc (scaled by its 



value at x pea k) using the scaling x/xhv introduced by Hurst & Vassilicos [15j |. One can clearly see that 
whilst use of xhv collapses the data obtained in the T=0.J^6m tunnel, a large discrepancy remains with 
the T=0.91m tunnel data. In particular, x pea k differs for BFG17 and SFG17. 
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Figure 4: (a) Schematic of wake interactions resulting from the fractal grid's bars, (b) Centreline 
mean streamwise velocity vs distance normalised by the wake- interaction length scale = L^/to 
(Uoo = 5.2m/s). 

The turbulence generated by either regular or multiscale/fractal grids with relatively low blockage 
ratio results from the interactions between the wakes of the different bars. In the case of fractal grids 
these bars have different sizes and, as a result, their wakes interact at different distances from the grid 
according to size andposition on the grid (see direct numerical simulations of turbulence generated by 
fractal grids in (l9| and (ri|). Assuming that the typical wake- width t at a streamwise distance 
x from a wake- generating bar of width/thickness tj (j 
such width corresponds to the largest bars on the grid, i.e. 



37|, the largest 



0, ...,N - 1) to be I ~ yjl 
i ~ \/tQX. Assuming also that this formula 
can be used even though the bars are surrounded by other bars of different sizes, then the furthermost 
interaction between wakes will be that of the wakes generated by the largest bars placed furthermost on 

~ ^t x±. 



the grid (see figure 4(a) ). This will happen at a streamwise distance x 
We therefore introduce 



x± such that L 



o 



x„ 



t2 

to 



(6) 



as a characteristic length-scale of interactions between the wakes of the grid's bars which might bound 
x peak- We stress that the assumptions used to define x+ are quite strong and care should be taken 
in extrapolating this presumed bound on x pea k to any space-filling fractal square grid beyond those 
studied here, let alone any fractal/multiscale grid. 

Figure [4(b)| is a plot of the normalised centreline mean velocity Uc/Uoo as a function of dimensionless 
distance x/x+ for all our four space- filling fractal square grids. All the data from the T- 0.46m tunnel 
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collapse in this representation. However the BFG17 data from the larger wind tunnel do not. They 
fall on a similar curve but at lower values of Uc/Uoq. This systematic difference can be explained 
by the fact that the air flow causes the BFG17 grid in the large wind tunnel to bulge out a bit and 
adopt a slightly curved but steady shape. The flow rate distribution through this grid must be slightly 
modified as a result. To compensate for this effect we introduce the mean velocity U p characterizing the 
constant mean velocity plateau in the vicinity of the fractal grids. In figure 5(a) we plot the normalised 
centreline mean velocity Uc/U p as a function of xjx* for all fractal grids and both tunnels and find an 
excellent collapse onto a single curve. 
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Figure 5: Streamwise evolution of (a) the centreline mean velocity Uq normalized by the initial mean 
velocity plateau U p and (b) centreline turbulence intensity normalised by its value at x = x pea k- Both 
plots are given as functions of x scaled by the wake-interation length-scale = L^/to (Uoc = 5.2m/s). 
In (b) the dashed line represents equation (9) with B = 2.06 and A = 2.82. 



In figure [5(b)] we plot u' c /Uc normalised by its value at x pea k as a function of x/x*. We also find an 
excellent collapse onto a single curve irrespective of fractal grid and tunnel. It is also clear from figures 
5(a) and 5(b)| that the longitudinal mean velocity gradient ^jr becomes insignificant where x > 0.6x* 
and that the streamwise turbulence intensity peaks at 



•Epeak 



0.45x, 



0.45^. 

to 



(7) 



The wake-interaction length-scale appears to be the appropriate length-scale characterising the first 
and second order statistics of turbulent flows generated by space-filling fractal square grids, at least on 
the centreline and for the range of grids tested in the present work. 



In figures 6(a) and 6(b) we plot the streamwise evolutions of the dimensionless centreline velocity 
Uc/Uoo and the centreline turbulence intensity u' c /Uc for various inlet velocities U^. These particular 
results have been obtained for the fractal grid SFG17 but they are representative of all our space-filling 
fractal square grids. One can clearly see that x pea k is independent of Uoq. Moreover, our data show 
that the entire streamwise profiles of both Uc/Uoq and u' c /Uc are also independent of the inlet velocity 
Uoq in the range studied. 



Hurst & Vassilicos |l5| have shown that the centreline turbulence intensity decays exponentially 
in the decay region x > x pea k and that the length-scale x pea k can be used to collapse this decay for 
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(a) 



x/M 



(b) 



x/M 



Figure 6: (a) Streamwise evolution of the centreline (y = z = 0) mean velocity in the T-O.46 tunnel 
for the fractal grid SFG17. (b) Streamwise evolution of the turbulence intensity for the same grid and 
in the same tunnel. 



different space- filling fractal square grids as follows: 



u'c _ u^{x ) ( , ( x-x 



ex P -B> (8) 



£/2 U%(x ) V V x peak 



where xo is a virtual origin and B' an empirical dimensionless parameter. We confirm this scaling decay 
form, specifically by fitting 
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uf / „ / x 



• I' sf; 



-4 = Aexp [-B - (9) 



to our experimental data with a slightly modified dimensional parameter B and an extra dimensionless 
parameter A which does not have much influence on the quality of the fit except for shifting it all up 
or down. We have arbitrarily set xq = 0, which we are allowed to do because the virtual origin does 
not affect the value of B. It only affects the value of A. 



As shown in figure 5(b) the exponential decay law (9) is in excellent agreement with our data for 
all the space-filling fractal square grids used in the present work. In particular, the parameter B seems 
to be the same for all the fractal square grids we used. Using a least-mean square method we find 
B = 2.06 and A = 2.82. 

Evidence in support of the idea that the decay region around the centreline downstream of Xp ea k is 



approximately homogeneous and isotropic was given in [3l| . The exponential turbulence decay observed 



in this region ([15l|. |3l| ) is therefore remarkable because it differs from the usual power-law decay of 
homogeneous isotropic turbulence. We have already reported in this subsection that the mean flow 
becomes approximately homogeneous in the streamwise direction where x > 0.6x*, i.e. in the decay 
region. In the following subsection we investigate the spanwise mean flow and turbulence fluctuation 
profiles downstream from space-filling fractal square grids and show how a highly inhomogeneous flow 
near the fractal grid morphs into a homogeneous one beyond 0.6x+. 



4.2 Homogeneity 

Like regular grids, the flow generated in the lee of space-filling fractal square grids is marked by strong 
inhomogeneities near the grid which smooth out further downstream under the action of turbulent 
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diffusion. This is evidenced in figures 7(a) and |7(bJ| which show scaled mean velocity U/U c and tur- 
bulence intensity u' /U profiles along the diagonal in the y — z plane, i.e. along the line parameterised 
by (y 2 + z 2 ) 1 / 2 in that plane. Close to the grid, the mean velocity profile is very irregular, especially 
downstream from the grid's bars where large mean velocity deficits are clearly present. These deficits 
are surrounded by high mean flow gradients where the intense turbulence levels reach local maxima as 
shown in figure [7(b)] Further downstream, both mean velocity and turbulence intensity profiles become 
much smoother supporting the view that the turbulence tends towards statistical homogeneity. Note 
that figures 7(a) | and |7(b)| show diagonal profiles at x/M e ff ~ 7 and 53 in the lee of the BFG17 grid 
which is a long way before i„ m t (see figure 3b). The profiles are quite uniform in the decay region as 
shown by Seoud Sz Vassilicos [3l| . Our evidence for homogeneity complements theirs in two ways: they 
concentrated only on the decay region whereas we report profiles in the production region and how 
they smooth out along the downstream direction; and we report diagonal profiles whereas the profiles 



m 



31| are all along the y coordinate. 
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Figure 7: (a) Normalised mean velocity profiles and (b ) turbulence intensity profiles along the diagonal 
line in the T-0.91m tunnel. The fractal grid is BFG17. Symbols: 0> X /M e ff ~ 7 (i.e. x/x* ~ 0.02); 
■ , x/M e ff « 53 (i.e. x/x* w 0.15). 



To evaluate the distance from the grid where the inlet inhomogeneities become negligible, we intro- 
duce the ratios U c /Ud and u' c /u' d where subscripts c and d denote respectively the centreline (y = z = 0) 
and the streamwise line cutting through the second iteration corner (y = z = 6cm in the SFG17 case) 



as shown in figure 8(a) These two straight lines meet the inlet conditions at two different points, their 
difference being representative of the actual inhomogeneity on the fractal grid itself: one point is at 
the central empty space whilst the other is at the corner of one of the second iteration squares. The 



streamwise evolutions of U c /Ud and u'JuL are reported in figure 8(b) As expected, large differences are 
observable in the vicinity of the space-filling fractal square grid for both the mean and the fluctuating 
velocities. One can see that the mean velocity ratio U c /Ud is bigger than unity. This reflects the mean 
velocity excess on the centreline where the behaviour is jet-like because of the local opening compared 
to the mean velocity deficit downstream from the second iteration corner where the behaviour is wake- 
like because of the local blockage. This difference between jet-like and wake-like local behaviours also 
explains why the fluctuating velocity ratio u' c /u' d is almost null close to the grid where the centreline 
is almost turbulence free. Further downstream, both ratios U c /Ud and u' c /u' d tend to unity as the flow 
homogenises. One can see that inhomogeneities become negligible by these criteria beyond x/x* = 0.6, 
which is quite close to x pea k- 

A main consequence of statistical homogeneity is that the small-scale turbulence is not sensitive to 
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Figure 8: (a) The measurements in (b) are taken in the T-0.J^6m tunnel along the straight lines which 
run in the x-direction and cut the plane of the pictured SFG17 grid at points c and d. (b) Streamwise 
evolution of the ratios U c /Ud (o) and u' c /u' d (*) measured for the SFG17 grid in the T-0.J^6m tunnel 
(Uoo = 5.2m/ s). The horizontal dashed lines represent the range ±10%. 



mean flow gradients. For this to be the case, the time scales defined by the mean flow gradients must 
be much larger than the largest time-scale of the small-scale turbulence. From our measurements and 
those in (3l| (see figure 3 in (3~D|). (dU / dx)~ l and (dU/dy) -1 are always larger than about 1 and 0.15 
second, respectively, at and beyond x pea k where the time-scale of the energy-containing eddies is well 
below 0.07 seconds. The ratios between the smallest possible estimate of a mean shear time scale and 
any turbulence fluctuation time scale are therefore well above 2 (in the worst of cases) at x pea k and far 
larger (by one or two orders of magnitude typically) beyond it for all three fractal grids SFG8, SFG13 
and SFG17 and all inlet velocities in the range tested. By this time-scale criterion, from x pea k 
onwards, the small-scale turbulence generated by our fractal grids, including the energy-containing 
eddies, is not affected by the typically small mean flow gradients which are therefore negligible in that 
sense. 

We close this subsection with figure 9 which illustrates in yet another way the homogeneity of the 
turbulence intensity at x > 0.5x* and suggests that equation (9) can be extended beyond the centreline 
in the y — z plane in that homogeneous region x > 0.5x*, i.e. 

£ = do) 

with the same values A ~ 2.82 and B ~ 2.06 independently of and space-filling fractal square 
grid. It is worth pointing out, however, that figure 9 also shows quite clearly that x pea k c an vary 
widely across the y — z plane and that a very protracted production region does not exist everywhere. 
Formula (7) gives x pea k on the centreline but x pea k can be very much smaller than 0.45x* at other 
y — z locations. This is a natural consequence of the inhomogeneous blockage of fractal grids and the 
resulting combination of wake-like and jet-like regions in the flows they generate. 



4.3 Skewness and flatness of the velocity fluctuations 



Previous wind tunnel investigations of turbulence generated by space-filling fractal square grids [15| . 



3l| have not reported results on the gaussianity/non-gaussianity of turbulent velocity fluctuations. We 



therefore study here this important aspect of the flow, mostly in terms of the skewness S u 
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Figure 9: Streamwise turbulence intensities as functions of x/x* along the centreline and along the 
straight streamwise line which crosses point d in Fig. 8(a). SFG17 grid in the T-0.J^6m tunnel with 
t/oo = 15m/ s. This plot remains essentially the same for the other values of Uqo that we tried. 

and flatness F u = ' 2 of the longitudinal fluctuating velocity component u along the centreline. This 
skewness is also a measure of one limited aspect of large-scale isotropy, namely mirror symmetry, as it 
vanishes when statistics are invariant to the transformation u to —u, but not otherwise. Isotropy was 



studied in much more detail in [15| and |31[ where x-wires were used. These previous works reported 
good small-scale isotro py in the decay region (3~0 | and levels of large-scale anisotropy before and after 
Xpeak on the centreline [15| which, for the turbulence generated by the grids SFG13, SFG17 and BFG17 
in particular, are very similar to the levels of large-scale anisotropy in turbulence generated by active 
grids 01, & 

We first check that both S u and F u do not depend on the inlet velocity Uqo, and this is indeed the 
case as shown in figures 10(a) and 10(b). These plots are particularly interesting in that they reveal the 
existence of large values of both S u and F u at about the same distance from the grid on the centreline. 
This distance scales with the wake-interaction length-scale x* as shown in figures 11 and 12. Indeed, 
the profiles of both S u and F u along the centreline collapse for all four fractal square grids (SFG8, 
SFG13 and SFG17 in the T-0.46m tunnel and BFG17 in the T-0.91m tunnel) when plotted against 
x/x*. The alternative plots against x/M e ff clearly do not collapse (see figures 11(a) and 12(a)). 

For comparison, figure 11(a) contains data of S u obtained with the regular grid SRG which are 
in fact in good agreement with usual values reported in the literature (see e.g. 0|, (26|). It is well 
known that regular grids generate small, yet non-zero, positive values of the velocity skewness S u and 



Maxey [2j| explained how their non-vanishing values are in fact a consequence of the free decay of 
homogeneous isotropic turbulence. Whilst the velocity skewness S u generated by fractal grids takes 
values which are also close to zero yet clearly positive in the decay region, S u behaves very differently 
in the production region on the centreline. 

The behaviours of S u and F u in the production region on the centreline are both highly unusual as 
can be seen in figures 10, 11 and 12. Clearly some very extreme/intense events occur at x ~ 0.2x* and it 
is noteworthy that the location of these extreme events scales with x* even though it is clearly different 
from x pea k ~ 0.45x* . The scatter observed at and around this location is due to lack of convergence 
because these intense events are quite rare as clearly seen in time traces such as the one given in figure 
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Figure 10: (a) Skewness and (b) flatness of the longitudinal fluctuating velocity component u as 
functions of x/M e ff along the centreline. SFG17 grid in the T-0.J^6m tunnel. The dashed lined in 
the right plot is F u = 3. 
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Figure 11: Velocity skewness S u as function of (a) x/M e jf and (b) x/x*. Uoo = 5.2m/ s. Data 
corresponding to all fractal grids and to the regular grid SRG as per insert. 



6, i 



\1 



(a) 



-O- SFG8 
-■- SFG13 
-V- SFG17 
-♦- BFG17 



V * Q 



100 

x/M „ 



o SFG8 

■ SFG13 

V SFG17 

« BFG17 



0.6 0.8 
X / X, 



Figure 12: Velocity flatness F u as function of (a) x/M e ff and f&j x/x*. [Too = 5.2m/s. Data corre- 
sponding to all fractal grids as per insert. The dashed line in both plots is F u = 3. 
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13(a). These intense events are so high in magnitude (\S U \ ~ 1) that they cannot be attributed to 
experimental uncertainties. The negative signs of S u and of these intense events on the time traces 
demonstrate clearly that these extreme events correspond to locally decelerating flow. We leave their 
detailed analysis for future study. 
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Figure 13: (a) and (b) Normalised velocity samples recorded downstream of the SFG17 fractal grid 
at (a) x ~ 0.2x* and (b)x ~ 0.72;*. (c) Probability Density Function (PDF) computed from signals 
corresponding to the SFG17 grid such as those shown on the left. The dashed line is the Gaussian 
distribution. Uoo = 5.2m/ s for all the plots on this figure. 

These extreme decelerating flow events are also reflected in the Probability Density Function (PDF) 
of u which is clearly non-gaussian and skewed to the left (i.e. towards negative u-values) at x = 0.2x* 
whereas it is very closely gaussian in the decay region (see figure 13(b)). The flatness F u takes values 
close to the usual gaussian value of 3 in the decay region and remains close to 3 for all x > x pea k (see 
figures 10(b) and 12(b)). 

Note, finally, that (10) and F u ~ 3 in the homogeneous region x > 0.5x* imply that 



<u 4 > 



in that region, with the same values A 
fractal square grid. 



3A 2 exp -2B 
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2.82 and B ~ 2.06 independently of and space- filling 



4.4 Spectral energy budget in the decay region 

In the region beyond x pea k where the turbulence is approximately homogeneous and isotropic, the 
energy spectrum has been reported in previous studies [3jJ, 15J to be broad with a clear power-law 
shaped intermediate range where the power-law exponent is not too far from -5/3. For this region we 
can follow George &; Wang [12j] who found a solution of the spectral energy equation 



0_ 

di 



E{k,t) = T(k,t) - 2vk 2 E{k,t) 



(12) 



that implies an exponential rather than power-law turbulence decay. In this spectral energy equation, 
E(k, t) is the energy spectrum and T(k, t) is the spectral energy transfer at time t. The energy spectrum, 
if integrated, gives |u /2 , i.e. 

o poo 

12 - / E(k,t)dk. (13) 
'o 



2 U 
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The correspondence between the time dependencies in these equations and the dependence on x in our 
experiments is made via Taylor's hypothesis. 



George & Wang [12j showed that (12) admits solutions of the form 

E(k,t) = E s (t)f(kl(t),*) (14) 

T(k,t)=T s (t)g(kl(t),*) (15) 

where the functions / and g are dimensionless, and that such solutions can yield an exponential turbu- 
lence decay. These solutions are single-length scale solutions (the length-scale l(t)) and therefore differ 
fundamentally from the usual Kolmogorov picture 0| which involves two different length scales, one 
outer and one inner, their ratio being an increasing function of Reynolds number. The argument * in 
(14) and (15) represents any dependencies that there might be on boundary/inlet/initial conditions. 

The exponential decay reported by George &; Wang |l2l | exists provided that the length-scale l(t) is 
independent of time, i.e. 4rl = 0, and takes the form u' 2 ~ exp [-10i/t/Z 2 ]. Unless/ 2 oc v, this form does 
not obviously compare well with the exponential decay (10) because (10) is independent of Reynolds 
number. We now carefully apply the approach of George & Wang (l2l | to our data by making explicit 
use of all potential degrees of freedom and confirm that an exponential decay which perfectly fits (10) 
can indeed follow from their approach. 

Consider 

E(k,t) = E s (t,U oo ,Re ,*)f(kl(t),Re ,*) (16) 

T(k,t) = T s (t,U oo ,Re ,*)g(kl(t),Re ,*) (17) 

where Reo = u ™ to is the Reynolds number which characterises the thickest bars on the fractal grid, 
l(t) = l(t, Reo,*) and the argument * includes various dimensionless ratios of bar lengths and bar 
thicknesses on the fractal grid, as appropriate. The functions / and g are again dimensionless. The 



conditions for (16) and (17) to solve (12) are (see [l2j for the solution method) 



d 2u 

-E, = -a^E, (18) 

T, = b f t E, (19) 

dl dE„ ,„ , . 

n' l = c nri E - (20) 

where a, b and c are dimensionless functions of Reo and * (note that a > 0). Under these solvability 
conditions, the spectral energy equation (12) collapses onto the dimensionless form 

d k 2 fin Reo *) 

f(n, Re , *) + c(Re , /(re, Re , *) = b(Re , *)g(n, Re , *) H { — (21) 

are a(Keo,*) 

where k = kl(t, Reo, *). 

Two different families of solutions exist according to whether c vanishes or not. If c 7^ 0, then c 
must be negative and 

l(t) = l(to)[l + ^(t-r )] 1 / 2 (22) 

E s (t) = E s (t )[l + ^(t - ro)] l ' 2c (23) 
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in terms of a virtual origin tq. However, if c = 0, then % = and 



It is this second set of solutions, the one corresponding to c = 0, with which George & Wang [12j] chose 
to explain the form (10). 

From (24), (16) and (13) and making use of t = x/Uoq it follows that 



u' 2 = u' 2 exp [~^J^ J (25) 

where u' = u'nQJ^, ReQ, *), a = a(ReQ,*) and I (Ren,*) is independent of x. Our wind tunnel 
measurements in subsection 4.1 and those leading to (10) and its range of validity suggest that (24) 
and (10) are the same provided that u'q/Uoo is a dimensionless function of the geometric inlet parameters 
* and nothing else, and that 

a = imRe Q l 2 (Re , *)/£§ (26) 



if use is made of (6) and B ~ 2.06. In other words, the single-scale solution of George & Wang |12| 
fits our data provided that the single length-scale I (ReQ,*) is independent of x (i.e. c = 0), that the 
dimensionless coefficient a in (18) depends on ReQ and * as per (26) and that u' scales with Uoq. Under 
these conditions, it follows from equation (21) that the dimensionless spectral functions / and g satisfy 

f(kl, Re Q , *) = b(Re , *)g(kl, Re , *) + 2R<% 1 ( -) (kl) 2 f(kl, Re , *) (27) 

\VBl(Re ,*)J 

where B ~ 2.06. Integrating this dimensionless balance over k = kl yields 



f°° ( l \ 2 r°° 

I f (k, Ren, *)dn = 2 ( —= I Re^ 1 / k 2 f (k, Ren, *)dn 

7o \VBl(Reo,*)J ° Jo 



(28) 



because the spectral energy transfer integrates to zero. This equality can be used in conjunction 
with (16), (13) and (6) to obtain a formula for the kinetic energy dissipation rate per unit mass, 
e = 2v J °° k 2 E(k,t)dk: 

QR„'2 rr n ,l2 T J 

J CO 



This is an important reference formula which we have been able to reach by applying the George & 
Wang [Til theory and by confronting it with new measurements which we obtained for three different 
yet comparable fractal grids and three different inlet velocities Uoo. These are the new measurements 
reported in subsections 4.1. and 4.2 and summarised by (10) and (6) along with the observation that A 
and B in (10) do not depend on [Too an d on the different parameters of the space- filling fractal square 
grids used. 

4.5 Multiscale-generated single-length scale turbulence 



No sufficiently well-documented boundary- free shear flow [34j , 29j nor wind tunnel turbulence generated 



by either regular or active grids [26l |. [27J] has turbulence properties comparable to those discussed here, 
namely an exponential turbulence decay (10), a dissipation rate e proportional to u' 2 rather than 
the usual u' 3 , and spectra which can be entirely collapsed with a single length-scale. It is therefore 
important to subject our data to further and more searching tests. 
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The downstream variation of the Reynolds number Re = is different for different boundary-free 
shear flows. However, it is always a power-law of the normalised streamwise distance x ^ x ° where Lb 
is a length characterising the inlet and xq is an effective/ virtual origin. For example, Re ~ ( x ^ x ° )" 1 ^ 3 
for axisymmetric wakes, Re ~ ( X 'lr" )° for plane wakes and axisymmetric jets and Re ~ ) 1 / 2 



for plane jets 34|, [29(|. The turbulence intensity's downstream dependence on 2_£o j s ^ ( x ^ x ° ) 2 / 3 



v Lb 

for axisymmetric wakes, ~ ( x ~[ Xo )" 1 / 2 for plane wakes and jets and ~ for axisymmetric 



jets 3J|, [29|. In wind tunnel turbulence generated by either regular or active grids the downstream 
turbulence also decays as a power law of and so does Re (26| . (2^ . In all these flows, as in 

fact in all well-documented boundary- free shear flows, the integral length scale L u and the Taylor 
microscale A grow with increasing x, and in fact do so as power laws of X L ~ X ° ■ Their ratio L u /X 
is a function of x ~ x ° and of an inlet Reynolds number Reg = Uoa ^ B where Uoo is the appropriate 



inlet velocity scale. Estimating A from vu /A ~ e ~ u /L u [34j . [29|, it follows that L u /X 



Rel /2 (^y 1/6 for axisymmetric wakes, L u /X ~ Rel /2 (^)° for pi ane wakes and axisymmetric 
jets, L u /X ~ Rel /2 (^-) l / i for plane jets and L u /A ~ ReJ 2 {^-)- p with p > 1/2 for wind tunnel 
turbulence. These downstream dependencies on Reo and x ^ x ° can be collapsed together as follows: 

^ ~ Rex (30) 

where Re\ = This means that values of L u /X obtained from measurements at different downstream 
locations x but with the same inlet velocity Uoo an d values of L u /X obtained for different values of 
Uoo but the same downstream location fall on a single straight line in a plot of L u /X versus Re\. This 
conclusion can in fact be reached for all sufficiently well-documented boundary- free shear flows (34| as 
for decaying homogeneous isotropic turbulence (2^ . 0| if the cornerstone assumption e ~ u' 3 /L u is 
used. 

The relation L u /X ~ Re\ is a direct expression of the Richardson-Kolmogorov cascade and, in 
particular, of the existence of an inner and an outer length-scale which are decoupled, thus permiting 
the range of all excited turbulence scales to grow with increasing Reynolds number. This relation is 
therefore in direct conflict with the one- length scale solution of George & Wang [l~2| . Seoud & Vassilicos 
(m found that L u /X is independent of Re\ in the decay region of the turbulent flows generated by our 
space-filling fractal square grids (in which case Reo is defined in terms of the wind tunnel inlet velocity 
Uoo an d Lb = to). Here we investigate further and refine this claim and also show that it is compatible 
with the theory of George Sz Wang [l~2l |. 

The essential ingredient in the previous subsection's considerations is the single-length scale form 
of the spectrum (16). Our hot wire anemometry can only access the ID longitudinal energy spectrum 
E u (k x ) of the longitudinal fluctuating velocity component u. The single-length scale form of E u (k x ) is 
E u (k x ) = E su f u (k x l, Reo, *) which can be rearranged as follows if use is made of ^u' 2 = J Q °° E u (k x )dk x : 

E u (k x ) = u' 2 lF u (k x l,Reo,*) (31) 

where F u (k x l, Reo,*) = \f u {k x l, Reo, *)/ J dk x lf u (k x l,Reo,*). In the case where u' 2 = u f2 (x,Uoo,*) 
decays exponentially (equations (24) to (29)), the length-scale / is independent of the streamwise dis- 
tance x from the grid. 

An important immediate consequence of the single length-scale form of the energy spectrum is that 



both the integral length-scale L u and the Taylor microscale A are proportional to I [12j, [3l| . Specifically, 



L u = l j dK x K x l F u (n x ,Reo,*) = al (32) 
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and 



A 




dK x K,lF u (K x ,Re , *) = (31 



(33) 



where a and (3 are dimensionless functions of Ren and *. This implies, in particular, that both L u and 
A should be independent of x (as was reported in 15 1 and (3l| ) if I is independent of x in the decay 
region. Using (29) and uu' 2 /X 



Uoo/x* we obtain 



A ~ LnRe, 



-1/2 



(34) 



which is fundamentally incompatible with the usual e 



~ u 
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/L u . As noted in [3l|, e ~ u' 3 /L u is in 



fact straightforwardly incompatible with an exponential turbulence decay such as (10) and an integral 
length-scale L u independent of x. 

We now report measurements of L u , A and E u {k x ) which we use to test the single-length scale hy- 
pothesis and its consequences. These measurements also provide some information on the dependencies 
of L,, and A on Ren and *. 



Firstly, we test the validity of (34). In figures 14(a) and 14(b) we plot X/\/ vx*/U versus x/x* along 



the centreline. These figures do not change significantly if we plot X/\/ ux^/Uoo = (A/ Ln)Re l J 2 versus 
x/x*. It is clear that (34) and scaling x with x* offers a good collapse between the different fractal 
grids where x > 0.2x*, and that A does indeed seem to be approximately independent of x/x* in the 
decay region as reported in (l5| and (3l| and as predicted by (33) and (34). However the collapse for 
different values of Uoo is not perfect and there seems to be a residual dependence on Ren which is not 
taken into account by (34). It is worth noting here that our centreline measurements for the regular 
grid SRG with Uoo = 5.2m/ s produced data which are very well fitted by (A/M e// ) 2 = 3 • 10~ 4 j;/M e// 
in agreement with previous results 0|, (26| and usual expectations [l|. 
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Figure 14: A/ \J 'vx*/U versus x/x* along the centreline, (a) Data for the four different fractal grids 
and Uoo = 5.2m/ s. (b) Data for the fractal grid SFG17 and three different inlet velocities Uoo- 

Before investigating Reynolds number corrections to (34), and therefore (29) which (34) is a direct 
expression of, we check that the turbulence generated by low-blockage space-filling fractal square grids 
is indeed fundamentally different from other turbulent flows. For this, we plot L u /X versus Re\ in 
figure 15(a) Whilst (30), which follows from e ~ u' 3 /L u , is very well satisfied in turbulent flows not 
generated by fractal square grids, it is clearly violated by an impressively wide margin in the decay 
region of turbulent flows generated by low-blockage space-filling fractal square grids. This is not just a 
matter of a correction to usual laws; it is a matter of dramatically different laws. 
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Figure 15: (a) L u /X as a function of R\ for different inlet velocities Uoo and different positions on 
the centreline in the decay region in the lee of the SFG17 fractal grid. The dashed line represents 
the empirical law obtained by fitting experimental data from various turbulent flows none of which 
is generated by fractal square grids. This is the data compiled in 25fl (which includes jet, regular 
grid, wake and chunk turbulence), and the dashed line is an excellent representation of L u /X ~ R\ 
particularly at Re\ > 150. (b) R\ as a function of x/x* for three different inlet velocities Uoo in the 
lee of the SFG17 fractal grid. 



One important aspect of (30) is that it collapses onto a single curve the x and ReQ dependencies 
of L u /X for many turbulent flows. Figures 16(a) and 16(b) show clearly that, in the decay region, 



L u /X is independent of x and also not significantly dependent on fractal square grid, but is clearly 
dependent on Uoo- There are other turbulent flows where L u /X is independent of x, notably plane 
wakes and axisymmetric jets. However, the important difference is that Re\ is also independent of 
x in plane wake and axisymmetric jet turbulence whereas it is strongly varying with x in turbulent 



flows generated by fractal square grids (see figure 15(b)). As a result, (30) holds for plane wakes and 
axisymmetric jets but not for turbulent flows generated by fractal square grids where, instead, L/X is 
independent of Re\ in the decay region (see figure 15(a) ), as previously reported in |3l| . It is not fully 
clear from figure 15(a) if L/X is or is not a constant independent of Rcq for large enough values of 
ReQ (specifically for values of Uoo larger or equal to 10m/ s in the case of 15(a)). The results in [3l| 



might suggest that L/X is independent of ReQ for large enough values of ReQ, but figure 16(b) does 
not comfortably support such a conclusion. More data are required for a conclusive assessment of this 
issue which is therefore left for future study. 

We now turn to the Reynolds number corrections which our measurements suggest may be needed 



in (34). Figures 17(a) and 17(b) show that 



A ~ LoRe, 



-1/3 



(35) 



is a better approximation than (34) as it collapses the different Uoo data better without altering the 
quality of the collapse between different fractal grids. From e ~ uu'^/X 2 , (35) implies 



u' 2 U 



e ~ 



"°Re~ 1/3 



(36) 



which is an important iZeo-deviation from (29) and carries with it the extraordinary implication that 
e tends to as ReQ — > oo. Of course, this implication is an extrapolation of our results and must be 
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Figure 16: L u /X as a function of x/x* along the centreline (a) for the four different fractal grids and 
Uoo = 5.2m/s and (b) for the fractal grid SFG17 and three different inlet velocities U^. 

dealt with care. In fact, we show in subsection 4.8, that this extrapolation is actually not supported 
by our data and the single-length scale theoretical framework of our work. 

There is only two ways in which George & Wang's fl2| single-scale theory can account for these 
deviations from (29) and (34). Either (i) these deviations are an artifact of the different large-scale 
anisotropy conditions for different values of Reo, or (ii) the single-length scale solution of (12) which 
in fact describes our fractal-generated turbulent flows belongs to the family for which c ^ 0, not the 
family for which c = (see subsection 4.4). 

(i) Large-scale anisotropy affects the dimensionless coefficient required to replace the scaling uu' 2 /X 2 r> 
e(= §-B - according to (29)) by an exact equality. This issue requires cross-wire measurements at 
many values of [Too to be settled and must be left for future study. 

(ii) If c / 0, i.e. c < 0, then 



U ' 



2E s {xq) \va\c 



\ [ 1 + 7ir^(---o)] (1 - c)/2c (37) 
(where we have used x = Uoot and xq = Uooto) and (32)-(33) remain valid but with 

l(x,Reo,*) = l(x ,Reo,*)[l + j^^ (x - x )] 1/2 (38) 



not with a length-scale I independent of x. Additionally, the following estimates for the Taylor mi 

d_ 

dt 



croscale and the dissipation rate e can be obtained, respectively, from ^u' 2 ~ uu' 2 /A 2 and from an 
integration over k of (21): 



A ~ l S^p4[l + j^-ix - x ^ 2 , (39) 
y/2a(l - c) l 2 {xo)Uoc y 

u' 2 

e = 3va(Re , *) (1 + \c(Re , *)|) (40) 

These two equations replace (34) and (29) which follow from c = 0. It is easily seen that the power-law 
form (37) tends to the exponential form (25) and that the Taylor microscale A becomes asymptotically 
independent of x — xq in the limit c — > 0. 
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The new forms (37)-(40) depend on two length-scales, 1(xq, Rcq, *) and Xq, one kinetic energy scale 
E s (xq)/1{xq) and two dimensionless numbers, va(Reo, *)/ (l(xo, Reo, *)£^x>) an d c(Reo,*), all of which 
may vary with Reo and boundary /inlet conditions. There seems to be enough curve- fitting freedom 
for these forms to account for our data in the decay regions of our fractal-generated turbulent flows, 
in particular figures 5b, 6b, 9, 14, 15, 16 and 17. In subsection 4.8 we present a procedure for fitting 
(37) and (39) to our data which is robust to much of this curve-fitting freedom. It is worth noting 
here, in anticipation of this subsection, that (39) is consistent with the observation (originally reported 
in [15l | and [3lJ) that A is approximately independent of x in the decay region, but provided that 
l^[xo)U ( x ~ x o) 1 m muc h of this region. However, (39) also offers a possibility to explain the 
departure from the constancy of A at large enough values of x/x* where A appears to grow again with 
x (see figures 14(a) and 17(a)), very much as (39) would qualitatively predict. Note, in particular, that 
this departure occurs at increasing values of x/x* for increasing Reo (see figures 14(b) and 17(b)) , 
something which can in principle also be accounted for by (39). In the following subsection we show 
that same observations can be made for L u . 
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Figure 17: A/ (Lo( ^*° )~ w ) versus x/x* along the centreline for n = 1/3 (various values of n were tried 
and n = 1/3 offers a clear best fit), (a) Data for the four different fractal grids and U^ = h.lmjs. (b) 
Data for the fractal grid SFG17 and three different inlet velocities U^. The results look very similar 
when plotting A/ (L (— )~ n )- 



4.6 The integral length-scale L u 

The streamwise evolution of the longitudinal integral length-scale L u on the centreline is plotted in 
figure 18 for all space-filling fractal square grids as well as for the regular grid SRG. The integral 
length-scales generated by the fractal square grids are much larger than the regular grid's even though 
their effective mesh sizes are smaller. Comparing data from the T-O.J^Gm tunnel, one can see that L u 
appears independent of the thickness ratio t r . However, there are large differences between the integral 
scales generated by fractal grids SFG17 and BFG11 which have the same t r but fit in different wind 
tunnels. This observation suggests that the large-scale structure of the fractal grid has a major influence 
on the integral length-scale. One might in fact expect that the integral length-scale is somehow linked 
to an interaction length-scale of the grid. For regular grids this interaction length-scale is typically the 
mesh size, whereas for space-filling fractal square grids a large variety of interaction length-scales exist, 



the largest being Lq. Figure 19(a) supports the view that the scalings of L u and its x-dependence are 



mostly determined by Lq and a;*, respectively, though not perfectly. Figure 19(b) suggests that L u /Lq 
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is not significantly dependent on the inlet Reynolds number Rcq, at least for the range of Uoo values 
investigated here. This figure was obtained for the SFG13 grid but is representative of our other three 
space-filling fractal square grids as well. 
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Figure 18: Centreline streamwise evolution of the longitudinal integral length-scale L u normalized by 
the mesh size M e ff (Uoo = 5.2m/s). Data obtained for the regular grid SRG are also shown (*). 
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Figure 19: L u /Lq as a function of x/x* along the centreline, (a) Four different fractal grids and 
Uoo = 5.2m/ s. (b) Three different inlet velocities Uoq with the SFG13 fractal grid. 

Irrespective of whether I does or does not depend on streamwise coordinate x, equations (32) and 
(33) suggest that L u /\ is definitely not a function of x but that it can nevertheless be, in all generality, 
a function of Reo and of the fractal grid's geometry. In fact figure 16 is evidence of some dependence 
on Uoq at least at the lower Uoq values. Assuming L u ~ Lq as seems to be suggested by figures 19(a) 

1 1 1 /2 1 

and 19(b) and using either (34) or (35) implies either L u /X ~ Re or L u /X ~ Re . Of these two 



implications, it is the latter which agrees best with our measurements (see figures 20(a) and 20(b) where 
we pi 
(34). 



1/3 

we plot (L u / A)/ Re versus x/x*) which is consistent with the fact that (35) fits our data better than 
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Figure 20: (L u /X) /Re as a function of x/x* along the centreline (a) for the four different fractal 
grids and Uoq = 5.2m/ s and (b) for the fractal grid SFG17 and three different inlet velocities [Too- 



As reported in [15l | and |31| . L u remains almost constant in the decay region (see figures 19(a) and 
19(b) ). The independence on x can be explained in terms of the c = single-scale solution of George & 
Wang (l2| . but it might also be even better explained in terms of the c < single-scale solution which 
yields (32) and (38), i.e. 

L u = a(Re , *)l(x , Re , *)[1 + ,^"^1 (x - x )] 1/2 . (41) 



A careful examination of figures 19(a) and 19(b) suggests that L u /Lq may be constant for some of the 
way downstream in the decay region till it starts increasing slightly, very much like the behaviour of 
the Taylor microscale A (see figures 14 and 17) and in qualitative agreement with (41). In fact, the 
ratio L u /X is predicted by (39) and (41) to be independent of x in the decay region, which is in full 
agreement with figures 16. 

The fact that L u /X is independent of x (figure 16) in the decay region where Re\ decreases rapidly 
with increasing x (figure 15b) is evidence against e ~ u' 3 /L u , not only in the context of c = single-scale 
solutions of (12), but also in the context of c < single-scale solutions of (12). Indeed, (37) and (41) 
are consistent with Uoo4^u' 2 ~ —u' s /L u only if c = —1, in which case Re\ = u'\/u is independent of x 



because of (37) and (39) and therefore in conflict with our experimental observations and those in |15| 



and [311 ] . In fact, the downstream decreasing nature of Re\ in the decay region of our fractal-generated 



turbulent flows imposes —1 < c < 0. 
4.7 The energy spectrum E u {k x ). 

The results of subsections 4.4 to 4.6 imply that the small-scale turbulence far downstream of low- 
blockage space-filling fractal square grids is either fundamentally different from the small-scale tur- 
bulence in documented boundary-free shear flows and decaying wind tunnel turbulence originating 
from a regular/active grid, or e ~ u' 3 /L u does not hold in these non-fractal-generated flows where the 
length-scale ratio L u /X is proportional to Re\ if one assumes e ~ u' s /L u . Our results therefore shed 
serious doubt on the universality of e ~ u' 3 /L Ul the cornerstone assumption present either explicitely 
or implicitely in most if not all turbulence models and theories [H], 0|, 29], 30], 



However, our data do not allow us to educe with full confidence a formula for the dissipation rate e 
in turbulence generated by space-filling fractal square grids. This issue is related to the fact that whilst 



26 



an exponential turbulence decay (equations 10 and 11) fits our data well, the Rcq dependence of A 
which follows from it in a self-preserving single- length scale context does not. Qualitative observations 
of the x-dependence of L u and A may suggest that the turbulence decay is in fact a power-law of the 
type (37), rather than exponential, albeit with a power-law exponent large enough (i.e. c close enough 
to 0) for the exponential form to be a good fit. An attempt at addressing this issue is made in the 
following and final subsection 4.8. This attempt relies on the results of our examination of energy 
spectra and the single-l eng th scale assumption which we now report. 

Seoud & Vassilicos [31[ studied the downstream evolution of the ID energy spectrum E u (k x ,x) in 
the decay region of space- filling fractal square grids and found that, for a given velocity Uoo, E u (k x ,x) 
can be collapsed for different downstream positions and for all our fractal grids in terms of (31) where 



/ is replaced by either A or L u . Indeed, E u {k x ) 



,'2 



\F u (k x X) and E u (k x 



L u F u (k x L u ) collapse 



the entire spectral data equally well at a given inlet velocity Uoo, a fact which we confirm in figures 
21(b), 22(a) and 22(b) These figures clearly support George Sz Wang's (l2| single-length scale 
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assumption (16) and (31). 
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Figure 21: SFG17 grid in the T-0.^6 tunnel with Uoo = 5.2m/s along the centreline. Compensated 
ID energy spectra normalized using (a) u' 2 and A or (b) u' 2 and L u . Spectra from different centreline 
positions x/x* collapse over all wavenumbers. 



However, Seoud & Vassilicos [3ll | did not attempt to collapse energy spectra for different inlet 
velocities Uoq. We therefore compare energy spectra obtained at the same position downstream of the 
same fractal grid but with different inlet velocities Uoq. In figures 23(a) and 23(b) we report such results 
obtained at x/x* ~ 0.62 with the SFG17 grid. These figures are representative of all other such results 
which we obtained with our other fractal grids and at other positions but which we do not present here 
for economy of space. 

From these figures, the form E u {k x ) 
different inlet velocities U^ than E u (k : 



u' 2 XF u (k x X) may seem to offer a much better collapse for 



a 



,2 L U F U { 



Z L U ). The discrepancy of the form E u ( 



,'2 



L u F u (k x L u ) is mostly at the high wavenumbers and is made evident by our compensation of the 
spectra by (k x L u ) 5 ^. It might be tempting to conclude that I is different from L u and in fact equal 
to A, but such a conclusion would be incompatible with E u (k x ) = u' 2 XF u (k x X) because (32) and (33) 
would then lead to the inconsistency that I is in fact not different from L u . The fact that L u /X grows 
with Reo, at least for the range of Rcq values considered here, suggests that we should be considering 

lF u (k x l, Re-o, *) without neglecting the dependencies on Rsq and perhaps 



a spectral form E u (k x 
even *. 

We have seen in the previous subsection that much of the scaling of L u is controlled by Lq, 



i.e. 
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Figure 22: T-O.46 tunnel with = 5.2m/s along the centreline. Compensated ID energy spectra 
measured at x/x* = 0.62 and normalized using (a) u' 2 and A or (b) u' 2 and L u . Spectra corresponding 
to our different grids appear to collapse over all wavenumbers. 



L u ~ Lq to a first scaling approximation, and that L u does not vary significantly with U^. Figure 
23(b) suggests that the plot of E u (k x )/ fu /2 L u ) versus k x L u is also imperceptibly dependent on Uqo at 



the lower values of k x L u but not at the higher ones. These three observations can all be explained if 
the assumption is made that 

l{x ,Re ,*) = L (42) 



and that 



F u (k x l,Re 



f u (k x l)H u {k x lRe^ 



(43) 



where n > and H u is a monotonically decreasing function which is very close to 1 where k x \R€^ n < 1 
and very close to where k x \Re^ n > 1. There may be residual dependencies on the geometry of the 
fractal grid, i.e. on *, but we do not have enough fractal grids in our disposal to determine them. Once 
again, this is an issue for future study. 

Equations (31), (38), (42) and (43) can readily account for the behaviour observed in figure [23(b)| 
Combined with (32), (42) and (43) also imply that L u scales with Lq provided that f u is a decreasing 
function of k x l where k x l > 1. Figures 20(a) and 20(b) would then suggest that /3 in (33) scales as 

1 /o 

Re . It is the function H u {k x lRe^ n ) in (43) which makes this scaling possible. In fact, if f u ~ {k x l)~ p 
where k x l > 1, then (33), (38) and (42) imply A ~ L Re^ n(3_p)/2 . Note that the Kolmo gorov-like 
exponents p = 5/3 and n = 3/4 (see [29| . (3-l |. 0|) yield A ~ L(jRe identically to (34) which follows 
from the c = single-length scale solution of the spectral energy equation (12). 



The good collapse in terms of both forms E u {k x ) 



u 
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XF u (k x X) and E u {k x ) = u u L u F u (k x L u ) in 



figures 21(a) and 21(b) comes from the fact that all data in these figures are obtained for the same 
value of Reo and the same fractal grid, and that L u /\ does not vary with x. These figures are therefore 
also consistent with (42) and (43). The good collapse of the form E u {k x ) = u' 2 L u F u (k x L u ) in figure 
22(b) is mainly a consequence of the fact that the fractal grids SFG8, SFG13 and SFG17 all have the 



same value of Lq and can also follow from (42) and (43). However, the apparently good collapse of 
the form E u (k x ) = u' 2 XF u (k x X) in figure 22(a) must be interpreted as being an artifact of the limited 

range of values of thicknesses tQ that we have experimented with (see table 1), more limited than the 

1 /3 

range of inlet velocities Uqo which allows the Re scaling of L u /X to be picked up by our spectra in 



figure 23(b) but not in figure 22(b) 
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Returning to figure 23(a) we notice that it does not, in fact, present such a good collapse of the 



data, particularly over the range of scales where the collapse in figure 23(b) appears good. Within 
the framework of (42) and (43), the semblance of a perhaps acceptable collapse in figure 23(a) results 
from a numerical circumstance to do with the exponents n and p. Chosing p = 5/3 and n = 3/4 for 
the sake of argument, (31), (38), (42) and (43) would imply that the quantity plotted in this figure, i.e. 
(k x \) 5 / 3 E u (k x ) I (u /2 A), is in fact equal to Re Q H u (k x \Re ) in the range which would correspond 
to k x X > 0.1 in the figure. Over the range of inlet velocities tried here, Re remains about constant 



— 1/3 

whilst Re varies a bit thus producing the effect seen in figure 23(a) a slight dependence on C/qo of 



the plateau and a semblance of a collapse of the dissipative range of the spectra. 
The conclusion of this data analysis is that the self-preserving spectral form 



E u (k x ) = u 2 lf u (k x l)H u (k x lRe ' 



(44) 



with 



r F1 4ua\c\ . 



xo)} 



1/2 



(45) 



is consistent with the theory of George Sz Wang [I2j and with our measurements in the decay region 
in the lee of our fractal square grids. We must stress again that future work is required with a wider 
range of fractal grids, measurement positions and inlet velocities in order to reach definitive conclusions 
confidently valid over a wider range of parameters. 
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Figure 23: SFG17 in the T-O.46 tunnel, (a) Compensated ID energy spectra normalized by the Taylor 
micro-scale A. (b) Compensated ID energy spectra normalized by the large-scale properties measured 
on centreline downstream fractal grid SFG17 for different inlet velocities Uoq {x/x* ~ 0.62 on the 
centreline) . 



4.8 Exponential versus power-law turbulence decay. 

We close section 4 with a discussion of how the exponential turbulence decay (25) corresponding to 
the c = solution and the power-law turbulence decay (37)- (38)- (39) corresponding to — 1 < c < 
solutions might fit together into a single framework. We already commented straight under equation 
(40) that the power-law form (37) tends to the exponential form (25) and that the Taylor microscale 
A becomes asymptotically independent of x — xq in the limit c — > 0. We now attempt to fit expression 
(37) to our decay region data and compare it with our exponential fit (10) (which is consistent with 
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(25) and (26)). To do this we start by fitting (39) to our Taylor microscale decay region data using 
some of the results reached in the previous subsection for our range of experimental parameters, namely 



/3 ~ Re ^ 3 and 1(xq) 



Lq. We therefore reformulate (39) as follows: 
A 2 



L 2 Re 



-2/3 



where 



7 



1 + 7 



4a |c| 
Re 



x 



X Q 



(46) 
(47) 

1/3 



We have, in effect, arbitrarily set to 1 a iieo-independent dimensionless parameter (or, if f3 ~ Re 
is somewhat faulty, perhaps even a weakly i?eo-dependent parameter) multiplying the left hand side of 
(46). However, most of the potential i?eo-dependence remains intact in this relation, in particular as 7 
is a priori i?eo-dependent. 

Relation (46) is plotted in Fig. 24(a) for the SFG17 grid and for all our inlet velocities [Too. From 
these curves it is easy to estimate 7 independently of the virtual origin xq, and we report our results 
in table 5. 



Uoo[m/s] 


5.2m/ s 


10m / s 


15m/s 


7 


0.23 


0.2 


0.13 



Table 5: Estimates of the coefficient 7 for the grid SFG17. 



We can now attempt to fit (37) to our turbulence decay data. Equation (37) can be recast in the 
form 



In (u' 2 ) = In 



2E s (x ) 
3l(x ) 



+ 



1 



2c 



In 



1 + 



Ava \ c\ 
l(x ) 2 U o 



(x - x ) 



(48) 



The observed near-constancies of A and L u in the decay region suggest from (39) and (41) that 

)re reaso 

2E s (xoY 



4l/a l c l (x — Xq) « 1. It is therefore reasonable to consider the first order approximation of (48) 



l(x ) 2 U, 
which is 



In (u' 2 ) w In 
and which can be reformulated as 



3Z(xo) 



+ 



1 



2c 



iua \ c\ 
l(x ) 2 U o 



(x - Xq) 



In (u' 2 ) w In 



2E s (x ) 
3l(x ) 



+ 5 



X — Xq 



(49) 



(50) 



with 5 = (-5-^) 7- This linear formula makes it easy to determine 5 from our experimental data 
independently of E s {xq) /1{xq) and xq, as indeed shown in figure 24(b) where (50) actually appears to 
fit our data well for all inlet velocities Uoq. Our resulting best estimates of the dimensionless parameter 
d are reported in table 6. This parameter appears to be iieo-independent, in agreement with the 



i?eo-independence of the turbulence intensity reported in figure 6(b) 



Uoolm/s] 


5.2m/s 


10m / s 


15m /s 


5 


-2.3 


-2.2 


-2.4 



Table 6: Estimates of the coefficient S for the grid SFG17. 
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The dimensionless coefficients a and c can now be obtained from our estimates of 7 and 5 using 

a = ^\5 + 1 /2\. (51) 

In table [7] we list the values thus obtained for c and a. It is rewarding to see that c turns out to be 
negative and in fact larger than — 1 . Of particular interest is the finding that c — > with increasing Reo 
and that the values of c are indeed quite close to for all our inlet velocities. These results suggest that 
the single-length scale power-law turbulence decay (37) tends towards the exponential turbulence decay 

(25) with the dimensionless coefficient a given by (51). Equation (51) is in fact equivalent to equation 

(26) which we obtained by fitting our turbulence intensity data with an exponential decay form. The 
Taylor microscale A also tends to an ^-independent form with increasing Reo because 7 — > 0, and so 
does 

L u = aL [l + 1 (x-x )/x^ 2 (52) 

(obtained from (41) and (42)). Indeed, we have checked that, in the decay regions of our fractal- 
generated turbulent flows, (52) provides a good fit of our L u data with the same values of 7 as the ones 
listed in table 5 and with a dimensionless constant (a ~ 0.34 in the case of the SFG17 grid) for all our 
inlet velocities U^. 

The dissipation rate e is given by (40) in the context of the power-law decaying single-length 
scale turbulence and it is easy to check that (40) tends to (29), the dissipation rate form of the 
c = exponentially decaying single-length scale turbulence, as Reo increases. Of course, this assumes 
that 7 and c tend to in that limit as the extrapolation of our fits would suggest. Equation (36) 
is incompatible with the view that power-law decaying single-length scale turbulence tends towards 

exponentially decaying single- length scale turbulence in the limit Reo — > 00. 

— 1/3 

Similarly, the empirical scaling of equation (35), i.e. A ~ LoRe , is also incompatible with such 
a gradual asymptotic behaviour. If use is made of (51), or equivalently (26), equation (39) shows 

— 1/2 

that, as Reo grows, A tends towards A ~ LoRe Q , the form predicted by the exponentially decaying 
single-length scale solution (see equation (34) and the argument leading to it). 

We noted in the previous subsection that an energy spectrum with a power-law intermediate range, 
i.e. f u ~ (k x l)~ p where k x l > 1, and a spectral form (31) with (43), (38) and (42) implies A ~ 

LoRe n ^ p ^ 2 . We also noted that the Kolmogorov-like exponents n = 3/4 and p = 5/3 yield A ~ 
—1/2 

LoRe . We are now suggesting that fits of the exponent n(3 — p)/2 might tend to 1/2 as Reo 
increases. This seems consistent with our observation that fits of the intermediate form f u ~ (k x l)~ p 
to our spectral data lead to p = 1.50 for Uoo = 5m/s, p = 1.57 for Uoo = lOm/s and p = 1.60 for 
Uoo = 15m/ s (see figure 25). The exponent p might indeed be tending towards 5/3 with increasing 
Reo, in which case we might also expect the exponent n to tend towards 3/4 if n(3 — p)/2 tends to 1/2. 



Uoolm/s) 


5.2m/s 


10m / s 


15m/s 


c 


-0.053 


-0.048 


-0.028 


a 


20.6 


24.6 


31.3 



Table 7: Estimates of the coefficients a and c for the grid SFG17. 
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Figure 24: SFG17 grid in the T-0.J^6m tunnel with three different inlet velocities, (a) Fit (46) of 
centreline A data, (b) Fit (50) of centreline u' 2 data. 



5 Conclusions and issues raised 

There are two regions in turbulent flows generated by the low blockage space-filling fractal square grids 
experimented with here. The production region between the grid and a distance about 0.5x* from it 
and the decay region beyond 0.5x*. In the production/decay region the centreline turbulence intensity 
increases/decreases in the downstream direction. The wake-interaction length-scale x* is determined 
by the large scale features of our fractal grids, = L^/to, but it must be kept in mind that one 
cannot change to an d or Lq without changing the rest of the fractal structure of these fractal grids. 
The downstream evolution of turbulence statistics scales on and can be collapsed with it for all our 
grids. However, it must be stressed that we have tested only four fractal grids from a rather restricted 
class of multiscale /fractal grids and we caution against careless extrapolations of the role of this wake- 



interaction length-scale to other fractal grids. For example, Hurst & Vassilicos [15j experimented with 
a low t r space-filling fractal square grid which seemed to produce two consecutive peaks of turbulence 
intensity instead of one downstream of it. A wider range of wake-interaction length scales should 
probably be taken into account for such a fractal grid, an issue which needs to be addressed in future 
work on fractal-generated turbulence. 

Whilst the turbulence in the production region is very inhomogeneous with non-gaussian fluctuating 
velocities, it becomes quite homogeneous with approximately gaussian fluctuating velocities in the decay 
region. Unlike turbulence decay in boundary-free shear flows and regular grid-generated wind tunnel 
turbulence where L u /X and Re\ change together so that their ratio remains constant, in the decay region 
of our fractal-generated turbulent flows L u /X remains constant and Re\ decreases as the turbulence 
decays. This very unusual behaviour implies that L u /X ~ Re\ and the Richardson-Kolmogorov cascade 
are not universal to all boundary- free weakly sheared/strained turbulence. In turn, this implies that 
e ~ u' 3 /L u is also not universally valid, not even in homogeneous turbulence as our fractal-generated 
turbulence is approximately homogeneous in the decay region. Inlet/boundary conditions seem to have 
an impact on the relation between L u /X and Reynolds number. The issue which is now raised for 
future studies is to determine what it is in the nature of inlet conditions and turbulence generation that 
controls the relation between the range of excited turbulence scales and the levels of turbulence kinetic 

1 /2 

energy. Whilst the general form L u /X = Reft fct(^^) may be universal, including fractal-generated 
turbulent flows, the actual function of X J^° in this form is not and can even be of a type which does not 
allow to collapse the x and Rcq dependencies by the Richardson-Kolmogorov cascade form L u /X ~ Re\. 

This issue certainly impacts on the very turbulence interscale transfer mechanisms, in particular 
vortex stretching and vortex compression which are considered to have qualitatively universal properties 
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Figure 25: (k x L u ) p E u (k x )/ (u' 2 L u ) versus k x L u for the value of the exponent p which gives the most 
horizontal plateau. Data obtained with the SFG17 grid at x/x* = 0.62 on the centreline of the T- 
0.46m tunnel with three different inlet velocities. The exponent p seems to increase towards 5/3 with 
increasing Uoo (p = 1.50 for U^ = 5.2m/s, p = 1.57 for U^ = 10m/ s, p = 1.60 for = 15m/s). 



such as the tear drop shape of the Q-R diagram [35j. Multi-hot wire anemometry [14j | applied to 
turbulence generated by low-blockage space-filling fractal square grids may have recently revealed very 
unusual Q-R diagrams without clear tear-drop shapes [171 ] . Fractal-generated turbulence presents an 
opportunity to understand these interscale transfer mechanisms because it offers ways to tamper with 
them. 

The decoupling between L u /X and Re\ can be explained in terms of a self-preserving single- 
length scale type of decaying homogeneous turbulence |l2l | but not in terms of the usual Richardson- 
Kolmogorov cascade ([l|, 0], (29| . (iol ]) and its cornerstone property, e ~ u' 3 /L u . This self-preserving 
single-scale type of turbulence allows for L u /X to increase with inlet Reynolds number Reo, as we in 
fact observe. This is a case where the range of excited turbulence scales depends on a global Reynolds 
number but not on the local Reynolds number. 

Our data support the view (both its assumptions and consequences) that decaying homogeneous 
turbulence in the decay region of some low-blockage space-filling fractal square grids is a self-preserving 
single-length scale type of decaying homogeneous turbulence [121 ] . Furthermore, our detailed analysis 
of our data suggests that such fractal-generated turbulence might be extrapolated to have the following 
specific properties at high enough inlet Reynolds numbers Re®: 

E u (k x ,x) = u ,2 (x)L (k x L )- 5 / 3 H u (k x L Re- 3/i ), (53) 

u' 2 {x) ^u%e- 2x/x % (54) 

e ss 3ii /2 [/oo/x*, (55) 

L u ~ L (56) 

and 

A ~ L Re~ 1/2 (57) 

where both L u and A are independent of x. A more detailed account of our conclusions involves the 
two types of single-scale solutions of the spectral energy equation, the c = and the — 1 < c < 
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types introduced in subsections 4.4 and 4.5. In subsection 4.8 we showed how our data indicate that 
the turbulence in the decay region is of the — 1 < c < type with a value of c which tends to as 
Re$ increases. This is why we stress the asymptotic extrapolations (53), (54), (55) and (57) in this 
conclusion. 

Our data require a very clear departure from the usual views concerning high Reynolds number 
turbulence 21 1, (29| . (3C)1 |. 34]. There is definitely a need to investigate these suggested high-i?eo 

properties further. Measurements with a wider range of fractal grids and a wider range of inlet velocities 
in perhaps a wider range of wind tunnels and with a wider range of measurement apparatus: x- wires, 



multi-hot wire anemometry [14| . [17| and particle image velocimetry. Direct Numerical Simulations 
(DNS) of fractal-generated turbulent flows are only now starting to appear [28| . fl~8| and their role will 
be crucial. Amongst other things, these studies will reveal dependencies on inlet/boundary geometrical 
conditions * which we have not been able to fully determine here because of the limited range of fractal 
grids at our disposal. 

A quick discussion of the features of extrapolations (53)- (57) reveals the various issues that they 
raise. The first issue which immediately arises is the meaning of Reo — > oo. We cannot expect this 
limit to lead to (53)- (57) if it is not taken by also increasing the number A of iterations on the fractal 
turbulence generator. How do our results and the extrapolated forms (53)-(57) depend on A? 

Secondly, in the extrapolated spectral form (53) we have assumed that the exponent p tends to 5/3 
in the high-i?eo limit and have therefore, in particular, neglected to consider any traditional small-scale 
intermittency corrections (see Q). This may be consistent with the observation of Stresing et al 33 1 



that small-scale intermittency is independent of Re\ in the decay region of our flows. However it is 
not clear why p should asymptotically equal 5/3 in the non-Kolmogorov context of our self- preserving 

—3/4 

single-scale decaying homogeneous turbulence. In particular, the inner length-scale LqRc differs 
from the Kolmogorov microscale (i^/e) 1 ^ 4 which scales as LqRcq 3//4 (io/^o) 1 ^ 2 (^'/^oo) _1//2 if account 

Q I A 

is taken of (55). If LoRe in (53) was to be replaced by this Kolmogorov microscale, then (57) 



would fail and the single-length scale framework of George & Wang [12j would fail with it. Why is the 
Kolmogorov microscale absent, or at least apparently absent, from decaying homogeneous turbulence 
in the decay region of some low-blockage space-filling fractal square grids? 

Thirdly, (55) suggests that the kinetic energy dissipation rate per unit mass is proportional to u' 2 
rather than u' 3 and that the turnover time scale is the global x*/Uoo rather than the local L u /u'. What 
interscale transfer mechanisms cause one or the other dependencies, and what are the implied changes 
in the vortex stretching and vortex compression mechanisms hinted at by the recent preliminary Q- 



R diagram results of Kholmyansky and Tsinober [17|? These issues directly address the universality 



questions raised in the Introduction and depend on the mechanisms of turbulence generation in the 
production region and the mechanisms which force important features of particular turbulence gener- 
ations to be or not to be remembered far downstream from the initial generator. What is the role of 
coherent structures, large or small, in shaping the type of homogeneous turbulence which decays freely 
in the decay region? 

Fourthly, is it possible that turbulence in various instances in industry and nature (e.g. in or over 
forest canopies, coral reefs, complex mountainous terrains, etc) might appear as a mixture of single-scale 
self-preserving turbulence and Richardson- Kolmogorov turbulence? Could such mixtures of two types 
of different turbulence give rise to what may appear as Reynolds number and intermittency corrections 
to the usual Richardson-Kolmogorov phenomenology and scalings? 

As a final note, it is worth comparing (55) with the usual estimate e = C e u' 3 /L u , which can also be 
seen as a general definition of the dissipation constant C e . One gets 

3 (to/Lp) 
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where use has been made of the estimate L u ~ 0.2Lq extracted from figures 19. The dissipation constant 
C t is not only clearly not universal, it can also be given bespoke values by designing the geometry of 
the turbulence-generating fractal grid, i.e. by changing the aspect ratio Furthermore, whilst 

a constant and universal value of C e would imply that, given a value of L u , the level of turbulence 
dissipation cannot come without an equivalent pre-determined level of turbulence fluctuations, (55) 
and (58) show that it actually is possible to generate an intense turbulence with reduced dissipation 
and even design the level of this dissipation. The implications for potential industrial flow applications 
are vast and include energy-efficient mixers (see Q) and lean premixed combustion gas turbines on 
which we will report elsewhere. 
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